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ABSTRACT 

The fraction of substructures required to account for anomalous flux ratios in gravita- 
tional lens systems appears to be higher than that predicted in the standard cold dark 
matter cosmology. We present a possible alternative route to anomalous flux ratios 
from lens galaxy environments. We consider the compound lens system such that a 
lens galaxy lie in a group or cluster, and estimate the contribution of substructures 
in the group/cluster to the fraction using an analytic model of substructures. We find 
that the contribution becomes dominant when the impact parameter of the lens is 
less than ^ 30% of the virial radius of the group/cluster. This indicates that the 
environmental effect can partly explain the high incidence of anomalous flux ratios. 
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1 INTRODUCTION 

It appears quite common that observed flux ratios between 
gravitationally lensed quasar images are difiicult to match 
while image posit ions can be fitt ed easily with simple smooth 
mass models. . Mao fc Schneideil ri998iil argued these anoma- 
lous flux ratios as evidence for substructures in lens galaxies. 
If this is true, it serves as a powerful probe of substruc- 
tures in galactic halos w h ich may not b e luminous enough 
("Mctcalf & Madau 2001; Chiba 2002; Dalai & Kochanck 
12002: .Keeton Gaudi. fc Fetters.. 200a: .Bradac et aL.2004) . 
For instance lPak^cKochmieM I2OO2I concluded that sub- 
structures should comprise 0.6% — 7% (90% confidence) of 
the mass of typical lens galaxies in order to explain flux 
ratios in several lens systems. Since the Cold Dark Mat- 
ter (CDM) model predicts 5% — 15% of mass in substruc- 
tures (e.g., iDe Lucia etal\ \2QQ4. IWeUer. Ostriker. fc Bod^ 
[200^, they argued that the result strongly supports for 
the CDM model. Not all anomalous flux ratios can be 
ascribed to the simplene s s of t he smooth pa r t of lens 
potentials (lEvans fc Witti 120031: iKawano etall \2004) or 
any propagation efi'ects because the degree of the (de- 
)magnification in observations depends on image parities 
fe.g.- iKochanek fc D alai 20oi) . 

However, a caveat is that the fraction of substruc- 
tures is in practice a very sensitive function of dis- 
tances froin__the_Jialo cente r: Both numerical simulations 
(" e.g.. iDe Lucia e t alJ l200'd) and analytic models (e.g.. 



lOeuri fc Lea 120041) clearly demonstrate that substructures 
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preferentially lie in outer parts of dark halos. This means 
that very little substructure is projected within the Einstein 
radius, since the Einstein radius is usually much smaller than 
the virial radius of halos. Thus the fraction of substructures 
in the CDM model appears to be sma^er than that required 
to explain anomalous flux ratios in gravitational lenses if 
we ta ke the spatial distribution of substr uctures into ac- 
count jEvans fc Witt'2003':' Mao et alJliopil . lMetcalj ll20ol . 
see also lchBi^Kravtsov. fc Keetonn2003F claimed that the 
discrepancy may be resolved by considering extragalactic 
halos along line of sight. Another possible explanation is 
some contribution s from stellar microlensing. For instance, 
IChiba et all ll2005h took mid-infrared images of B1422-I-231 
and FG1115-I-080, and concluded that the anomalous flux 
ratio of FG1115-I-080 is likely to be caused by stellar mi- 
crolensing. 

In this Letter, we present a possible alternative route 
to anomalous flux ratios. Specifically, we consider the 
compound lens system such that a lens galaxy lies in a 
group or cluster, and discuss to what extent substructures 
in the group/cluster could contribute to cause anoma- 
lous flux ratios. Although the situation is not universal, 
recent studies indicate that such compound systems 

(including incidenc e of foreground groups) appe ar to 

be quite common iKeeton. Christlein. fc ZabludoflJ l200d: 
i Blaridford. Surpi. fc Kundid I2OOII: iKeeton fc Zabludofj 
12004 iFassnacht et alJ 120041: lOeuri et al.l |2005^ . Obser- 
vations imply that many of them lie in groups, and 
lens systems in clusters are less common. The efl^ect of 
correlated matter outside the l ens halo was studied by 
IChen. Kravtsov. &: KeetonI i2003l) , however they considered 
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typical field galaxies as lens objects, thus the situation is 
somewhat difi'erent from ours. Throughout the Letter we 
assume a Lambda-dominated cosmology with the mass 
density parameter fim = 0.3, the vacuum energy density 
parameter Q.a = 0.7, the dimensionless Hubble constant of 
h — 0.7, and the rms fluctuation normalization as — 0.9. 
The results presented here will be, however, depend little 
on the specific choice of cosmological parameters. 



MASS FRACTION OF SUBSTRUCTURES 
FROM AN ANALYTIC MODEL 



Following iMao et alj ll2004ll . we compute the fraction of the 
projected surface density, /sub, to discuss the efficiency of 
flux anomalies. To do so, we use an analytic model of sub- 
structures developed by lOeuri fc Led ll2004) . The model 
takes account of two dominant dynamical processes that 
drive dominantly subhalo evolution: one is mass-loss caused 
by tidal interaction with the host halo, and the other is 
the orbital decay caused by dynamical friction which drives 
massive subhalos to the inner part of the host halo. The 
model also considers the formation epoch variation of the 
host halo, and the orbital decay of satellite halos outside 
the host halo virial radius. It was shown that the de- 
rived analytic distributions agree well wit h the results of 
recent high-resolution N-body simulation dPe Lucia et alJ 
12004). Although the model may be less realistic than semi- 
analytic approache s 'Bullock. Kravtsov. & Weinberg 2000; 
i Havashi etaP EoO.1: .Tavlor fc Babul .2004: .Zentner et aL 
l2005l: iPenarrubia fc Bensonll2005ft . it has the advantage of 
being able to compute spatial and mass distributions very 
easily. 

The fraction of the projected surface density fsuh is de- 
fined by the projected mass in substructures within an an- 
nulus with radii R- and R+ divided by the projected total 
mass within the annulus. We choose R- and R+ so that 
the annuli are equally spaced in log(ii/rvir) with a step size 
0.2, where r^ir is a virial radius of a halo that corresponds 
to a lens galaxy. First, we compute /sub in an isolated halo 
with mass M. The mass inside a cylinder with radius R are 
computed as 



47r 



M{< R) = ^ rf(cos 9)p{< r„,ax; M)r 



rmax = mm 



(r- 

V smb' 



(1) 
(2) 



where p(< r;AI) is the average density within a sphere of 
radius r, and 6 is an angle with respect to the projection di- 
rection. We can calculate the density profile of substructures 
Psub and corresponding mass Msub(< R) by integr ating the 
distributions of substructures llQguri fc Led 120041) over the 
mass of substructures. Then the mass fraction of substruc- 
tures is 



/sub — 



AMsub(fi-,i^+) 
AM{R^,R+) 



(3) 



where the mass within annuli is denoted as AM{R- , R+) = 
M(< R+)-M{< R-). 

Next we consider the situation that a lens galaxy lies 
in a group or cluster. In this case, the lens galaxy should 
correspond to a substructure (hereafter referred as primary 
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Figure 1. Calculation of /sub for a compound lens. The pro- 
jected masses M(< R), A/h(< R), and Afjj sub(< R) are given by 
equations 0, Jl}, and jSJ, respectively. 

substructure) in a host halo with mass Mh . Then to compute 
the lens potential we have to consider the mass associated 
the host halo (group or cluster) as well as the mass of the 
substructure. Here we consider the situation that another 
substructure in the host halo causes anomalous flux ratios 
in gravitational lenses (hereafter referred as secondary sub- 
structure). To quantify this contribution, we compute /sub 
in this compound lens system. When R is much smaller than 
the virial radius of the host halo, the mass of dark matter 
and substructures associated with the host halo can be ap- 
proximated as 

Mh(<i?) ^ 2^i?2 / '"■%(r;A'4)^=L=dr, (4) 
Jb V - 

Mh,sub(<i?) ~ 2-nR' / psub(r;Mh)^==dr, (5) 

Jb yf - &2 

with rh,vir being the virial radius of the host halo with mass 
A/h and h being the impact parameter of the primary sub- 
structure (see Figure 0. Then the mass fraction from sub- 
structures in the host halo is 



/sub — 



AMh,sub(J?-,J?+) 



l\M{R-,R+) + AMi,{R-,Rj, 



(6) 



where AM(_R_ , R+) is the mass associated with the primary 
substructure calculated from equation Q . Although the pri- 
mary substructure suffers from tidal stripping which makes 
the size of the substructure smaller, we neglect this because 
the contribution of the outer part to the mass estimate is 
quite small. 



3 RESULT 

We calculate /sub in the following setup: The mass of the iso- 
lated halo or the primary substructure, which corresponds 
to the lens object, is M = 5 x \Q^^h~^ Mq. As a redshift, we 
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Figure 2. The mass fraction of substructures /sub as a func- 
tion of projected radius R. The fractions of an isolated lens (eq. 
PI) and compound lens (eq. |^) are plotted by solid and dashed 
lines, respectively. We assume the lens (sub-)halo with M = 
5 X 10^2 /i-^M©, the host halo with = lO^'^h'^MQ, the im- 
pact parameter of the lens subhalo in the host halo b/rji^vir = 0-2. 
They are placed at 2 = 0.3. The vertical dotted lines i ndicate the 
posit ion of typical lensed images, 

■R/j-vir = 1% - 3% <Mao et alj 



choose a typical lens redshift of z = 0.3. For the compound 
system, we set the mass of the host halo Mh = lQ^'^h~^ Mq 
(i.e., Mh/M — 20). The virial radii of the lens and host halos 
are rvir ~ 295ft~^kpc and rh.vir ~ 800ft~^kpc, respectively. 
The primary substructure is placed so that the impact pa- 
rameter becomes fo/r-h,vir = 0.2. The value of the impact pa- 
rameter is slightly smaller than the average impact param- 
eter calculated from the radial distribution of substructures 
assuming the ratio of substructure to ho st halo masses to be 
10"^ (see Figure 3 of lOeuri fc Le3l2004l . 6/rh,vir ~ 0.3, but 
is a reasonable value if we take account of the bias due to 
the convergence of the host halo (see We integrate the 
mass of substructures in the range 10"" < m/M < 10"^ 
(5 X W^h'^MQ < m < 5 X 10" / i~^Mp ) to compare 
our results with those of iMao et alJ l)2004h . However, al- 
most uni versal form of the mass function of substructures 
Te3l2004h assures that relative amount of /a. 



m (V: L/eeli:^UU4ll assures that relative amount ot /aub be- 
tween isolated and compound lens systems does not change 
very much even if we shift the range to lower masses. 

Figure |5| shows /aub as a function of the projected ra- 
dius R. We plot /aub for both isolated compound systems. 
We note that for the compound system we show /sub from 
host group/cluster only: The lens galaxy in group/cluster 
itself should also have substructures that will increase /sub 
of the compound lens system further. For the isolated halo, 
we reproduce the result of JAap et al^ ( 20o3l . We find that 
/aub for the compound lens is comparable to, or even higher 
than, that for the isolated lens at the position of typical 
lensed images, R/r^ir ~ 1% — 3%. To examine how the re- 




Figure 3. Dependences of /sub for the compound lens on several 
parameters: The impact parameter b (upper left), the redshift 
z {lower left), and the mass of the host halo Mb {lower right). 
We change each parameter by fixing the other parameters to the 
fiducial values. 



suit is dependent on our specific setup describe above, in 
Figure 13 we show the dependence of /sub for the compound 
lens on several parameters. We find that /sub is quite sen- 
sitive to the impact parameter b but rather insensitive to 
other parameters including the mass of the host halo. There- 
fore we conclude that substructures in the host halo of the 
compound lens could be a dominant source of anomalous 
flux ratios when the impact parameter is 6/rh,vir 0.3. The 
predicted values, however, are still /sub ^ lO"'^ (without 
baryons) and thus they do not completely eliminate, though 
do reduce, the discrepancy between predicted and observed 
fraction of substructures. 

The expected level of the contribution from substruc- 
tures in a group/cluster can also be estimated by a simple 
argument. In our fiducial parameter set, the convergence at 
image position {R/vvir ~ 1% — 3%) is k ~ 0.25 (assuming 
a source redshift Zs = 1), and the external convergence is 
Kenv ~ 0.05. On the other hand, at 20% of a virial radius 
the fraction of substructures is ~ 0.06 (see Figure I^J. From 
these numbers, /aub is estimated as 

0.06 X Kcnv 



/aub(i?/rvir= l%-3%) 



' 0.01. 



(7) 



The estimated number is roughly consistent with the result 
presented in Figure |5| The external convergence Konv is a 
strong function of the impact parameter, thus the reason 
that /aub is sensitive to the impact parameter b can be also 
understood from this expression. 
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4 SUMMARY AND DISCUSSION 

We have examined the impact of lens galaxy environments 
on substructure lensing. Using an analytic model of sub- 
structures constructed bv lOeuri fc Leel l)2n04f) , we have com- 
puted the fraction of substructures in the situation that a 
lens galaxy lies in a group or cluster. Specifically, we have 
regarded a subhalo in a halo as a lens galaxy, and esti- 
mated the contribution of other substructures in the halo 
to the fraction of substructures which determines the effi- 
ciency of substructure lensing. We have found that the frac- 
tion depends mainly on the impact parameter of the pri- 
mary substructure (lens galaxy); the contribution becomes 
dominant when the impact parameter is less than ~ 30% 
of the virial radius of the host halo. We emphasize that 
the effect is rather insensitive to the mass of the host halo; 
our result suggests that even a poor group could have a 
great impact on substructure lensing. Therefore we con- 
clude that lens galaxy environments are important when 
we discuss anomalous fiux ratios. This, in turn, means that 
lens systems in dense environments may be more desir- 
able sites for detecting dark halo substructures, with ei- 
ther anomalous flux ratios or more direct methods (e . g . , 
lYonehara. Umemura. fc Susall2003l : llnoue fc Chi'ball2005ll . 

The result appears to c ontradict with that of 
IChen. Kravtsov. fc KeetorJ i2003r) who showed that nearby 
correlated structure is highly sub-dominant compared to 
substructure within the lens halo. However, they assumed 
average galaxies for lens objects and estimated the number 
density of clumps outside the lens using two-point corre- 
lation function ^(r), n ~ n{l + ^(r)), where n is the av- 
erage number density of clumps in the universe. In con- 
trast, we have assumed rather specific environment such that 
lens objects correspond to galaxies in a group or cluster, 
thus the situation is somewhat different from theirs. In our 
case, it is more appropriate to estimate the number density 
with the number density of clumps inside the group/cluster, 
fiG ~ Avir??- (where Avir ~ 300 is the characteristic value 
of the overdensity of collapsed objects), instead of n. This 
means that the number density of clumps in our situation 
should be higher by a factor ~ Avir, and this extra factor 
accounts for the difference between our and their results. 

The importance of this new route ma y be enhanced fur - 
ther if we consider the effect of baryon. As lMao et alJ (l2004l) 
noted, the baryon infall makes the density profile of lens 
galaxies steeper, p oc r~^, and hence substructures will be 
disrupted more easily. However, the argument is not applied 
to compound lens systems because the impact of baryon 
cooling on groups and clusters is much less significant. 

It is i mportant to address how ubiquitous s uch lens 
systems are. iKeeton. Christlein. fc Zabludofll i2000ll predict 
that ^ 25% of lens galaxies are i n groups or clusters, and 
iBlandford. Surpi. fc Kundid i200 1*) argued that the fraction 
coul d be even higher . Alth ough iPremadi fc M artel ( 2 0o3) 
and iDalal fc WatsonI i2005ft found somewhat smaller envi- 
ronmental effect, it appears that the models predict signifi- 
cantly smaller values for external shear than those obtained 
by modeling strong lensing data. In addition, it should be 
taken into account that anomalous flux ratios have been 
discussed primarily using quadruple lens systems. This in- 
dicates the sample would be biased significantly toward 
the higher fraction of the compound systems, because en- 



Table 1. Summary of known environments of lens systems 
with anomalous flux rat i os. T he sa mple is constru cted from 
iKeeton. Gaudi. fc Patter j |20o3) and iKochanek fc D alai |20oi). 
Note that the anomalous flux ratios of some of these system may 
be caused by stellar microlensing rather than lensing by substruc- 
tures. 



Name 


Environment 




MG0414-f0534 






B0712-I-472 


erouD (Fassnacht & 


Lubin 2002) 


RX J0911-I-0551 


cluster fBade et al. 


1997) 


SDSS J0924+0219 






PC 1115+080 


erouD fKundic et al 


. 1997) 


RX J1131-1231 






B1422-I-231 


erouD fTonrv 1998") 




B1933-(-503 






B2045-(-265 


erouD fFassnacht et 


al. 2004) 



vironmental effects enhanc e the quadruple to double ratio 
iKeeton fc Zabludo9l2004f) . Another effect we have to think 
of is the bias associated with the convergence of the group 
or cluster which increases lensing cross sections by a factor 
~ (1 — k)~^. These effects preferentially select compound 
systems with small impact parameters, indicating that the 
compound lens system considered in this paper may be quite 
ubiquitous. Indeed, the survey of environments of lens sys- 
tems with strong fiux anomalies, which is summarized in 
Table shows that more than half of the lens systems lie in 
groups/clusters, implying a correlation between flux anoma- 
lies and dense galaxy environments. 
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